Resonantly enhanced tunneling of Bose-Einstein condensates in periodic potentials 
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We report on measurements of resonantly enhanced tunneling of Bose-Einstein condensates loaded 
into an optical lattice. By controlling the initial conditions of our system we were able to observe 
resonant tunneling in the ground and the first two excited states of the lattice wells. We also 
investigated the effect of the intrinsic nonlinearity of the condensate on the tunneling resonances. 
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Resonantly enhanced tunneling (RET) is a quantum 
effect in which the probability for tunneling of a particle 
between two potential wells is increased when the quan- 
tized energies of the initial and final states of the process 
coincide. In spite of the fundamental nature of this ef- 
fect [i[ and the practical interest [2] , it has been difficult 
to observe experimentally in solid state structures. Since 
the 1970s, much progress has been made in construct- 
ing solid state systems such as superlattices [1, |3, Q and 
quantum wells [Q| which enable the controlled observa- 
tion of RET 0. 

In recent years, ultra-cold atoms in optical lattices Q 
have been increasingly used to simulate solid state sys- 
tems. Optical lattices are easy to realize in the labora- 
tory, and their parameters can be perfectly controlled 
both statically and dynamically. Also, more compli- 
cated potentials can be realized by adding further lattice 
beams [9| . This makes them attractive as model systems 
for crystal lattices, and in the last few years cold atoms 
and Bose-Einstein condensates (BECs) in optical lattices 
have been used to simulate phenomena such as Bloch os- 
cillations [l0[ and the Mott insulator transition [ll[. In 
this Letter we show that BECs in accelerated optical lat- 
tice potentials are ideally suited to studying RET. While 
in solid state measurements of RET only a few poten- 
tial wells were used and the periodic structures had to 
be grown for each realization, in our experiment the con- 
densate is distributed over several tens of wells and the 
parameters of the lattice can be freely chosen. Moreover, 
we are able to control the initial conditions of the system 
and thus observe RET in any chosen energy level and can 
also add nonlinearity to the system. 

A schematic representation of RET is shown in Fig. 1 . 
In a tilted periodic potential, atoms can escape by tun- 
neling to the continuum via higher-lying levels. The tilt 
of the potential is proportional to the force F acting on 
the atoms, and in general the tunneling rate T^z can be 
calculated using the Landau- Zener formula How- 
ever, when the tilt-induced energy difference FdL^i be- 
tween wells i and z -I- A? matches the separation between 
two quantized energy levels, the tunneling probability is 
resonantly enhanced and the Landau-Zcner formula no 
longer gives the correct result, as previously investigated 
in [13[ for cold atoms in optical lattices. While for the 
parameters of our experiment the enhancement over the 
Landau-Zener prediction was around a factor of 2 (see 
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FIG. 1: Explanation of resonantly enhanced tunneling. Tun- 
neling of atoms out of a tilted lattice is resonantly enhanced 
when the energy difference between lattice wells matches the 
distance between the energy levels in the wells. 



theoretical and experimental results of Fig. [Ha)), in gen- 
eral it can be several orders of magnitude. 

The starting point of our experiments is a BEC of 
^''Rb atoms, held in an optical dipole trap whose frequen- 
cies can be adjusted to realize a cigar-shaped condensate. 
The BECs are created using a hybrid approach in which 
evaporative cooling is initially effected in a magnetic 
time-orbiting potential (TOP) trap and subsequently in 
a crossed dipole trap. The dipole trap is realized using 
two intersecting gaussian laser beams at 1030 nm wave- 
length and a power of around 1 W per beam focused to 
waists of 50 /im. After obtaining pure condensates of 
around 5 x 10* atoms the powers of the trap beams are 
adjusted in order to obtain an elongated condensate with 
the desired trap frequencies (« 20 Hz in the longitudinal 
direction and 80 — 250 Hz radially) . 

Subsequently, the BECs held in the dipole trap are 
loaded into an optical lattice created by two gaussian 
laser beams (A = 852 nm) with 120 /im waist inter- 
secting at an angle 9. The resulting periodic poten- 
tial V{x) — Vb sin^(7ra;/dL) has a lattice spacing — 
X / {2 sm{9 / 2)) and its depth Vq is measured in units of 
the recoil energy E^cc = ?i^7r^/(2md|), where m is the 
mass of the Rb atoms. In the present experiment, we 
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used dL = 0.426 /zm (for ^0^00 = 6, 4, 9 and 16) and 
dL = O.Q20ijm (for Vq/E,^^ = 2.5, 10, 12 and 14). By 
introducing a frequency difference Aj/ between the two 
lattice beams (using acousto-optic modulators which also 
control the power of the beams), the optical lattice can 
be moved at a velocity v = dL^v or accelerated with an 
acceleration a = dL^^- 

A ramp from to Vq in around 1 ms loads the BEC 
adiabatically into the optical lattice [l3|. For loading 
the ground state levels, the lattice velocity is w = dur- 
ing the ramp. For the first and second excited levels, 
during the ramp the lattice is moved at a finite velocity 
calculated from the conservation of energy and quasimo- 
mentum [l^. Finally, the optical lattice is accelerated 
with acceleration a for an integer number of Bloch oscil- 
lation cycles. In the rest frame of the lattice, this results 
in a force F — ma on the condensate. Atoms that are 
dragged along by the accelerated lattice acquire a larger 
final velocity than those that have undergone tunneling, 
and are spatially separated from the latter by releasing 
the BEC from the dipole trap and lattice at the end of 
the acceleration period and allowing it to fall under grav- 
ity for 5 — 20 ms. After the time-of- flight, the atoms are 
detected by absorptive imaging on a CCD camera using 
a resonant flash. 

From the dragged fraction N drag /N tot, we then deter- 
mine the tunneling rate r„ in the asymptotic decay law 

Ndragit) = Ntot exp {-T,,t) (1) 

where the subscript n indicates the dependence of the 
tunnnelling rate on the local energy level n in which the 
atoms are initially prepared (ground state: n = I, first 
excited state: n = 2, etc.). In the experiments reported 
in this work, the number of bound states in the wells was 
small (2-4, depending on the lattice depth), so after the 
first tunneling event, the probability for tunneling to the 
next bound state or the continuum was close to unity. 

The resolution of our tunneling measurement is given 
by the minimum number of atoms that we can distin- 
guish from the background noise in our CCD images, 
which varies between 500 and 1000 atoms, depending on 
the width of the observed region. With our condensate 
number, and taking into account the minimum accelera- 
tion time limited by the need to spatially separate the two 
fractions after time-of-flight and the maximum accelera- 
tion time limited by the field of view of the CCD camera, 
this results in a maximum Tn/v-cec of « 1 and a minimum 
of « 1 X 10^^, with the recoil frequency v^cc — E^cc/h. 

A typical plot of the tunneling rate Fi out of the 
ground state as a function of Fq^ (where Fq = EdL/E^cc 
is the dimensionless force) in the linear regime is shown 
in Fig. 2(a). This regime is reached either by choos- 
ing small radial dipole trap frequencies or by releasing 
the BEC from the trap before the acceleration phase and 
thus letting it expand. In both cases, the density and 
hence the interaction energy of the BEC is reduced. Su- 
perimposed on the overall exponential decay of Ti/Fq 
with Fq^, one clearly sees the resonant tunneling peaks 
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FIG. 2: Tunneling resonances in an accelerated optical lattice, 
(a) Tunneling resonances of the n = 1 lowest energy level for 
Vo = 2.5 Erec- The arrows indicate the upper and lower limits 
for our precise measurement of r„. Inset: Deviation from the 
Landau-Zener prediction. For clarity, in both graphs only 
one representative error bar is shown, (b) Positions of the 
Ai = 1 resonance peaks as a function of the lattice depth. 
Only data points for which the resonance is clearly visible 
(e.g. not Ai = 1 of (a)) are included. Inset: Positions of the 
peaks for Ai = 2 and 3. 

corresponding to Ai — 2,3 and 4 (for this choice of pa- 
rameters, the Ai — 1 peak lay outside our experimen- 
tal resolution). In order to highlight the deviation from 
the Landau-Zener prediction, in the inset of Fig. [2] (a) 
we plot Ti/Tlz, where the Landau-Zener tunnelling rate 
Tlz is given by [Hill 

TLZ^I^recFoe . (2) 

The experimental results are in good agreement with nu- 
merical solutions obtained by diagonalizing the Hamilto- 
nian of the open decaying system [13, Il8l | Figure ^h) 
summarizes our results for the positions of the ground- 
state resonances Ai = 1,2 and 3 as a function of the 
lattice depth together with a theoretical fit assuming the 
separation of the lowest energy levels to be 

AE = aE,,,^/Vo/E,,,. (3) 
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Independently of Ai, the best fit is achieved for a = 1.5, 
to be compared with a — 2 for the harmonic osciUator 
approximation. A value a < 2 is to be expected since 
our lattice wells only contain a few bound states and are, 
therefore, highly anharmonic. 




FIG. 3: Anti-crossing scenario of the RET rates, (a) Theo- 
retical plot of ri,2 for Vb = 2.5-Erec with experimental points 
for Fi. (b) Theoretical plot of Fi,2 for Vb = 10i?rec with ex- 
perimental points for F2. For clarity, the vertical axes have 
been split and the r„ plotted on a linear scale, and only one 
representative error bar is shown. 

Using BECs in optical lattices allows us to explore res- 
onant tunneling in regimes that are difficult or even im- 
possible to access in solid state systems. First, we can 
prepare the condensates in the excited levels of the lat- 
tice wells before the acceleration. Again, tunneling res- 
onances are clearly visible, and the experimental results 
agree with theoretical calculations. The accessibility of 
higher energy levels allows us to experimentally deter- 
mine the decay rates at resonance of two strongly cou- 
pled levels. Although our experimental resolution does 
not allow us to measure the decay rates in two differ- 
ent levels for the same set of parameters Fq and Vq, we 
are able to compare the ground and excited state decay 
rates Fi and F2 with the theoretical predictions for two 
different parameter sets, as shown in Fig. 3. This figure 
reveals the anti-crossing of the decay rates of strongly 
coupled levels as a function of our control parameter Fq. 
These results demonstrate a peculiar behaviour of the 
Wannier-Stark states studied theoretically [gl, [l^ and 
more recently rephrased within a general context of cross- 
ings and anti-crossings for the real and imaginary parts 
of the eigenvalues of a non-hermitian Hamiltonians [20j . 
Our data confirm the predictions of 0] that the anti- 
crossings modify the decay rates of the two perturbing 
states in different ways. 

Additionally, by exploiting the intrinsic nonlinearity 
of the condensate due to atom-atom interactions, we can 
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FIG. 4: Resonant tunneling in the nonlinear regime, (a) Res- 
onance Ai = 3 for Vb ~ 2.5-Erec with C — 0.024 (squares), 
C = 0.035 (circles) and C = 0.057 (triangles). The solid line 
is the theoretical prediction for C = 0; the dashed lines are 
guides to the eye. (b) Dependence on C of the tunneling rate 
at the position of the peak Fq~^ = 1.21 (solid symbols) and 
of the trough _Fg~^ = 1.03 (open symbols). The dashed lines 
are fits to guide the eye. For clarity, in (a) and (b) only one 
typical error bar is shown. 



study RET in the nonlinear regime, as simulated in [2l|. 
In order to realize this regime, we carry out the accel- 
eration experiments in radially tighter traps (radial fre- 
quency > 100 Hz) and hence at larger condensate densi- 
ties. Figure 4(a) shows the results for increasing values 
of the nonlinear parameter (2^ 



(4) 



where uq is the peak condensate density and the s- 
wave scattering length. Two effects are visible: First, 
the overall (off-resonant) level of Ti increases linearly 
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with C. This is in agreement with our earher exper- 
iments on nonlinear Landau-Zener tunnehng p^ . [23| 
and can be explained describing the condensate evolu- 
tion within a nonlinearity-dependent effective potential 
Kff = Vq/{1 + 4C) [11]. Second, with increasing nonhn- 
earity, the contrast of the RET peak is decreased and the 
peak eventually vanishes. This is confirmed by the differ- 
ent dependence on C of the on- and off-resonant values 
of Fi (Fig. 4 (b)). We estimate that in order to signifi- 
cantly affect the resonant tunneling rate, the nonlinearity 
parameter has to be comparable to the width of the RET 
peak. This order-of-magnitudc argument agrees with our 
observations. 

Finally, we have experimentally tested the robustness 
of RET against a dephasing of the lattice wells induced by 
non-adiabatic loadin g of the BEC into the lattice in the 
nonlinear regime [l^, 25 1. Even for completely dephased 



wells, the tunneling resonances survive. 

In summary, we have measured resonantly enhanced 
tunneling of BECs in accelerated periodic potentials in a 
regime where the standard Landau-Zener description is 
not valid. Our results in the linear regime agree with nu- 
merical calculations, and the possibility to observe RET 
for arbitrary initial conditions and parameters of the pe- 



riodic potential underlines the advantage of our system 
over solid state realizations. Furthermore, we have ex- 
plored RET in the nonlinear regime and demonstrated 
that, as theoretically predicted, the tunneling resonances 
disappear for large values of the nonlinearity. 

In the present set-up the measurement of the tunnel- 
ing rate is limited in its dynamic range by the detection 
geometry. A larger dynamic range can be realized by 
long-distance transport of BECs 26]. Our method for 
observing RET can also be generalized in order to study 
other regular or disordered potentials, the effects of noise 
and the presence of a thermal fraction in the condensate. 
Furthermore, one might exploit the tunneling resonances 
to explore the spatial decohercncc processes and to per- 
form precision measurements. 
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